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Abstract 

In this paper we prove factorization of fragmentation function in non-equilibrium QCD by using 
Schwinger-Keldysh closed-time path integral formalism. We use the background field method of 
QCD in a pure gauge in path integral approach to prove factorization of fragmentation function in 
non-equilibrium QCD. Our proof is valid in any arbitrary gauge fixing parameter a. This may be 
relevant to study hadron production from quark-gluon plasma at high energy heavy-ion colliders 
at RHIC and LHC. 
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I. INTRODUCTION 



RHIC and LHC heavy-ion colliders are the best facilities to study quark-gluon plasma 
in the laboratory. Since two nuclei travel almost at speed of light, the QCD matter formed 
at RHIC and LHC may be in non-equilibrium. In order to make meaningful comparison 
of the theory with the experimental data on hadron production, it may be necessary to 
study nonequilibrium-nonperturbative QCD at RHIC and LHC. This, however, is a difficult 
problem. 

Non-equilibrium quantum field theory can be studied by using Schwinger-Keldysh closed- 
time path (CTP) formalism However, implementing CTP in non-equilibrium at 
RHIC and LHC is a very difficult problem, especially due to the presence of gluons in 
non-equilibrium and hadronization etc. Recently, one-loop resummed gluon propagator in 
non-equilibrium in covariant gauge is derived in |3|, |4 1 . 

High pt hadron production at high energy e^e~, ep and pp colliders is studied by using 
CoUins-Soper fragmentation function js-?]. For a high px parton fragmenting to hadron, 
Collins-Soper derived an expression for the fragmentation function based on field theory and 
factorization properties in QCD at high energy. This fragmentation function is universal 
in the sense that, once its value is determined from one experiment it explains the data at 
other experiments jsj. 

Recently we have derived parton to hadron fragmentation function in non-equilibrium 
QCD by using Schwinger-Keldysh closed-time path integral formalism |91]. This can be 
relevant at RHIC and LHC heavy-ion colliders to study hadron production from quark- 
gluon plasma. We have considered a high px parton in medium at initial time tq with 
arbitrary non-equilibrium (non-isotropic) distribution function f{p) fragmenting to hadron. 
The special case f{p) = — corresponds to the finite temperature QCD in equilibrium. 

e T ±1 

We have found the following definition of the parton to hadron fragmentation function 
in non-equilibrium QCD by using closed-time path integral formalism. For a quark (g) with 
arbitrary non-equilibrium distribution function fq{k) at initial time, the quark to hadron 
fragmentation function is given by js] 
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where z (=-^) is the longitudinal momentum fraction of the hadron with respect to the 
parton and Pt is the transverse momentum of the hadron. In the above equations \%n> is 
the initial state of the non-equilibrium quark (gluon) medium. The path ordered exponential 

= V exp[z^ /° d\n- A\x^ + n^\) T*^] (2) 

J —oo 

is the Wilson line 

Eq. ([1]) can be compared with the following definition of Collins-Soper fragmentation 
function in vacuum 
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Eq. ([T]) differs from eq. (I3]) only by the presence of prefactor 
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and by the replacement 



of the vacuum |0 > by the state \in >. 

Similar to the pp collisions, the evolution of this non-equilibrium fragmentation func- 



ll|-ll3| which 



tion at RHIC and LHC can be studied by using DGLAP evolution equation 
we have recently extended to non-equilibrium QCD jl^ . 

Another important physics analysis one needs to do at RHIC and LHC is to prove factor- 
ization of fragmentation function in non-equilibrium QCD to study hadron production from 
quark-gluon plasma. Factorization refers to separation of short distance from long distance 
effects in field theory. 

In this paper we will prove factorization of fragmentation function in non-equilibrium 
QCD. Factorization of fragmentation function in pp collisions is studied by Collins, Soper 
and Sterman by using diagrammatic technique to all orders in perturbation theory. Recently, 
we have proved factorization of soft and coUinear divergences by using background field 
method of QC D |l5| in a pure gauge in path integral approach. We have used the gauge 
fixing identity [l6| which relates the generating functional Z[A, J,r],f]] in the background 
field method of QCD in pure gauge with Z[J,ri,f]] in QCD (without the background field). 

We find it convenient to use background field method in a pure gauge to prove factoriza- 
tion of soft and collinear divergences in non-equilibrium quantum field theory [l7|. In this 



paper we will extend this to non-equilibrium QCD. Unlike QED |10|, the gauge fixing term 



and ghost term in the background field method of QCD [l8|, ll9| depend on the background 
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field Hence the extension of proof of factorization from non-equilibrium QED to 

ibrium QCD is not straightforward. For this reason we will use the gauge fixing 
16[ to prove factorization of fragmentation function in non-equilibrium QCD by 



non-equi 
identity 



using Schwinger-Keldysh closed-time path integral formalism in background field method of 
QCD in a pure gauge. Our proof is valid in any arbitrary gauge fixing parameter a. This 
can be relevant to study hadron production from quark-gluon plasma at RHIC and LHC. 

The paper is organized as follows. In section II we briefly review the factorization of 
soft and collinear divergences in QCD in pp collisions by using background field method 
of QCD. In section III we describe Schwinger-Keldysh closed-time path integral formalism 
in non-equilibrium QCD relevant for our purpose. In section IV we prove factorization of 
fragmentation function in non-equilibrium QCD by using background field method of QCD 
and closed-time path integral formalism. Section V contains conclusions. 



II. BACKGROUND FIELD METHOD OF QCD AND FACTORIZATION OF 
SOFT AND COLLINEAR DIVERGENCES 



Background field method of QCD was originally formulated by 't Hooft [18| and later 
extended by Abbott [19]. This is an elegant formalism which can be useful to construct gauge 
invariant (off-shell) green's functions in QCD. This is because, unlike QCD, the lagrangian 
density in the generating functional in the background field method of QCD is background 
gauge invariant (even after quantizing the theory). In the presence of external sources a 
relation between the generating functional Z[J, r^,?/] in QCD and the generating functional 



Z\A^ J, ?7, r/] in the background field method of QCD in pure gauge is recently obtained 16 |. 
We will use this relation to prove factorization of fragmentation function in non-equilibrium 
QCD. 

In QCD, the generating functional is given by 

^QCD[J,r7,r7] =/[rfQ][#][^iV^] det(^^^) e^/^^-[-i^"^.[Q]-2^(9.Q^")^+'^^[Q]'^+^-Q+^'A+.v;]. 

(4) 

Under the infinitesimal gauge transformation the quantum gluon field transforms as 

^Ql = -gr'^u^'Ql + d,u^ (5) 
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where is the gauge transformation parameter. In the background field method of QCD 



the generating functional is given by 
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^background qcd[A J , V , 7]] = I [dQ][dip][dip] det(- 
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i J d^x[-^F'^l,[A+Q]-^{G^{Q)f+^pmA+Q]i,+J-Q+fnP+v^] ^g^ 



where is the background field. The gauge fixing term is given by 

G'^(Q) = d.Q^'' + gr'^A'^Q^' = D^[A]Q'^'^ (7) 
which depends on the background field A" The infinitesimal gauge transformation is given 



by (l9| 



5AI = 0. (8) 
Changing the variable Q — !■ Q — ^ in eq. (EI) we find 

^background qcb[A, J, f], f]] = e"^/'^'--^-^ J [dQ] [#] [#] det(^^) 

/ d^xl-lF-l^lQ]-^{G'}{Q))^+'4>mQ]i'+J-Q+v^+vi^] (9) 

where 

G%Q) = d^Q^'' + gr'^A'^Q'^' - d^A'^'^ = D^[A]Q^'' - (10) 

The most crucial statement of factorization theorem of Collins, Soper and Sterman is the 
appearance of the Wilson line eq. (|2]) in the definition of the parton distribution function 

El 

1 ^...-^'"V„„ 



^troirac ^tr,oior[7+ < ^lV^(y",Z/T) , Vt] $"^0] ^(0)|P >]. (11) 

and in fragmentation function in eq. ([3]) which makes them gauge invariant. This Wilson 
line is responsible for cancelation of soft and collinear divergences which arise due to presence 
of loops and/or higher order Feynman diagrams. This Wilson line can be thought of as a 
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quark or gluon jet propagating in a soft and/or collinear gluon cloud. The Wilson line eq. 
([2]) can be written as 

= V exp[ig r dXh- e^^-^A^ix") T"] = V exp[ig-^h ■ A%x^) T^]. (12) 

J— oo h • o 

Since the vector is free we can choose it to correspond to various physical situations. For 
example if we choose = n^, where is a fixed lightlike vector, we find from eq. ( !T2|) the 
phase factor 

uj'^(x) = -L^n-A\x). (13) 
n ■ a 

Using the Fourier transformation 

^U^) = J (04^M(^)e^''^ (14) 



V = g uj{k) = ig (15) 



we find from eq. (IT3 



Note that eq. ( !T5|) is precisely the Eikonal vertex for a soft gluon with momentum k inter- 
acting with a high energy quark (or gluon) jet moving along the direction js, 6|. In the 
soft gluon approximation in |6|] the fixed lightlike vector is taken to be having only "+" or 
"-" component: 

= (n+,n-,nT) = (1,0,0) or = (n+, n", nr) = (0, 1, 0). (16) 

Now we will show that when the classical background field is an abelian-like pure gauge 

Alix) = d.co^ix) (17) 

we can reproduce eqs. (fT3|l and (fT5|) which appear in the Wilson line eq. (fT2l) . Multiplying 
x'^ independent four vector h^^ from left in eq. f[T7|) we find 

h-A%x) = h-duj''{x). (18) 

Dividing h ■ d from left in the above equation we find 

co^(x) = j^h-A^ix). (19) 
h ■ o 
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Since the four vector is free we can choose it such way that A^^ can represent Feynman 
rules involving soft gluons or coUinear gluons. For example, when /i^ = n'^ where n^^ is a 
fixed light like vector (see eq. (ITB]1 ). we find from the above equation 

uj'^^x) = ^n-A^ix) (20) 

which reproduces eq. f|T3|) which appears inside the Wilson line in eq. (1121) . 
Similarly, if we choose = n^, where is a non-light like vector 

< = K,^B,0), (21) 

we reproduce the Feynman rules for the collinear divergences 6|]. This establishes the 
correspondence between the Wilson line and the classical background field A^^{x) as given 
by eq. (ITTj) in the context of soft and/or collinear divergences. 

The non-perturbative correlation function in the presence of the background field A'^^{x) 
is given by 

S S 

T— 7 r T 7 r ^background QCd[-4, J, 7] , f]]\ j=^=fj=o = < 1pi{x2)lpj{xi) >background QCD (22) 

Similarly the non-perturbative correlation function in QCD (without the background field) 
is given by 

5 5 

= <^.(a;,)^,(x,)>^5^ (23) 

When the background field A^(a;) is SU(3) pure gauge in QCD given by 

T''A''^ = -{d^U)U-\ u^^igT-^-(x) (24) 

'^9 



we find from 



15| 



< iJi{x2)i^j{Xi) >Lkground QCD = [P f d\ h ■ ^'^(x^ + X)T'']] 

J —oo 

_ _ P — OO 

X [< A{x2Wxi) >^'^'b] X [Vexp[zg dX h ■ A'ix-^ + h-" X)T% (25) 



All the Al^{x) dependences (responsible for soft and collinear divergences) have been factored 
into the path ordered exponentials or Wilson lines. This proves factorization of soft and 
collinear divergences by using background field method of QCD. 

Note that the structure function and fragmentation function are proportional to the 
non-perturbative correlation function in eq. ( l25|l with the following identification. In the 
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generating functional Z[A, J,r],f]] in the background field method of QCD the integration 
is only over those functions of quantum gluon field Q^(a;), quark field ipi{x) and antiquark 
ield ipi{x) such that their Fourier transform ipi{k) and ipi{k) vanish in the soft region 



l3- 



III. CLOSED-TIME PATH INTEGRAL FORMALISM IN NON-EQUILIBRIUM 
QCD 

Unlike pp collisions, the ground state at RHIC and LHC heavy-ion collisions (due to the 
presence of a QCD medium at initial time t = tm (say tj„=0) is not a vacuum state |0 > 
any more. We denote \in > as the initial state of the non-equilibrium QCD medium at tin- 
Consider the time evolution of the density matrix 

i^ = [Hj,p], p(-oo)=po (26) 

where Hj is the interaction hamiltonian. The formal solution is 

p(t) = S{t, -oo)po5(-oo, t), S{t, -oo) = T exp[~i f dt' Hi{t')]. (27) 

J —oo 

The density matrix p is in interaction picture. The average value of an operator L in the 
interaction picture is given by 

< Lit) >= Tr[pit)L{t)], = [Ho, Lit)] (28) 

where Hq is the free hamiltonian. Since in many situations we deal with correlation function 
of several fields at different times, it is useful to transfer all the time dependence to field 
operators and consider the density operator as independent of time, i.e. to go to Heisenberg 
representation. For the time independence of the density matrix we can take the value 
of the matrix determined by expression eq. f l27|) at a certain fixed instant of time, for 
example, t = tm = 0, having thus included in it all the changes which the distribution po 
had undergone when the external field and the interaction in the system were switched on. 

Consider the medium average of T— products of several operators. Using the Heisenberg 
density matrix one finds 

< T[L(t)M(t')-.] >= Tr[pT[Lit)Mit')...]] = Tr[S'(0, -oo)po5'(-oo, 0)T [L it) M it' )....]] 
= Tr[poSi-oo, 0)T[L(t)M(t')--]5(0, -oo)]. (29) 



Going over to the operators in the interaction picture 

< T[L{t)M{t')...] >= Tr[poS{-oo, 0)T[5(0, t)Li{t)S{t, t')Mi{t')....]S{Q, -oo)] 
^Tr[pon[SMtWiit')....]] (30) 

where Tc is the complete contour from 

- oo ^ i ^ t'.... -oo (31) 

and Sc is the complete yS— matrix defined along Tc. 

To deal with Feynman diagram and Wick theorem it is useful to split the time interval 
to "+" and "-" contour; where "+" time branch is from — oo to +00 where (time) T— order 
product apply and "-" time branch is from +00 to —00 where (anti-time)T— order product 
apply 

Consider scalar field theory first. Since there are two time branches there are two fields 
and two sources and hence four Green's functions. Let us denote the field (f)'^{x) and the 
source J'^{x) in the "+" time branch and (f)'{x) and J~{x) in the "-" time branch. The 
generating functional is given by 

Z[p, J+, J_] = / [#+][#-] exp[i[5[0+] - + / d'^xJ+(l)+ - [ d^xJ_0_]] < 0+, 0|p|0, 0_ 



where S[(l)\ is the full action in scaler field theory and |0±, > is the quantum state corre- 
sponding to the field configuration (f)±{x, t — 0). 

In the CTP formalism in non-equilibrium there are four Green's functions 

G^^{x,x') = .ffJt'wf!^ inmx)cf>{x')\in >=< r0(x)0(x') > 

G^4x,x') = , S^/r/l iri\fcf>{x)cl>{x')\in >=< fcf>{x)cf>{x') > 

[—Ij'^d J -[X) J -[X ) 

G+_{x,x') = _^^^j^(^J^j =< in\(l){x')(j){x)\in >=< (j){x')(l){x) > 

G_+{x, x') = tf/'^/fV^^ =< in\(Pix)(P{x') \in >=< 0(x)0(x') > (33) 

— I OJ—yXjJ-^-yX J 

where T is the time order product and T is the anti-time order product given by 

T^{x)4){x') = e{t - t')(t){x)(t){x') + e{t' - t)4){x')4){x) 

f(t){x)<t){x') = e{t' - t)(t){x)(t){x') + e{t - t')(t){x')(t){x). (34) 
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A. Generating Functional in Non-Equilibrium QCD 



Additional ghost fields are present in the gauge theory. However, we will directly work 
with the determinant of the gauge fixing terms. The generating functional in non-equilibrium 
QCD is given by 



exp[z / d\[-\{F^l^[Q^] - F'^lAQA) - ^{{d,QT? - {d.Q'-f) + 

X < g+, ^/'+, o| p |o, i)-,i).,Q-> (35) 

where p is the initial density of state. The state ■i/'^, t/'^, > corresponds to the field 
configurations Q'^{x,t = tin = 0), ■ip{x,t = tin = 0) and ip{x,t = tin = 0) respectively. Note 
that we work in the frozen ghost formalism j^, 4| for the medium part at the initial time 

IV. PROOF OF FACTORIZATION IN NON-EQUILIBRIUM QCD 

The generating functional in the background field method of QCD is given by eq. ([9]) 
where the gauge fixing term Gf{Q) is given by eq. ( ITO|) . Now combining the background 
field method of QCD with the closed-time path integral formalism we find the generating 
functional in non-equilibrium QCD in the presence of background field 

^background QCD [p, A, J+, J_, 7]+, f]+, r]_] = g"^ / '^*-(-^+-^+--^-^-) 

60"" (O ) SC^iO ) 

X det( ^ ) X det( /V"^ )x 

exp[^ J d'x[-\{F-l[QA - F'^llQ^]) - ^{{{G}{Q^)f - {G){Q-)f) + 

-V'-i/fg-]^- + J+-Q+- j--Q- + n+-iJ+-n--^- + v+-i'+-v-- i'A] 

X < Q+, 0| p |0, Q- > . (36) 

The fermion fields and corresponding sources transform as follows 

= U_ = ij_ UZ\ ri_ = U_ r/_, C = V- UZ^ (37) 
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where 

U+ = e^^^^^-T^-), U. = e'sT'^'^-^x) _ (3g) 

When the background field is pure gauge in QCD given by eq. fl24|) we find from fl5| 

ZQCD[J,r/,r/] =e*/'^'--^-^ x ^background qcd[A ^, r/', r/'] - j [dQ][di,][d^] det{^-^j^) 

exp[^ j d'x[--^F'^lM - ^iG%Q)Y + VI/[g]^ + J ■ g + # + v^]] X [t J d'xJ ■ SQ + ...] 

(39) 



where we have used the gauge fixing identity 16| . Since we work in the frozen ghost formal- 
ism ^, 4| for the medium part at the initial time t = tin = 0, the initial state > 
associated with the medium part at initial time is gauge invariant by construction. Hence 
we find by extending eq. (139|) to non-equilibrium QCD 



ZqcdIp, J+, V+, V-,V+, V-] = e*/'^'^[-^+-^+ -^'-^-^ X ^background qcd[p, A, J+, J_, r/;, r]'_, r]'^, r]'_] 
- j[dQ^][dQAW^mA[di^^]WA X det( £^+0 X det(^^) x 
x[i j d^x [J+ ■ 6Q+ - J_ ■ + ...] X 

exp[^ / d'x[-\{F-l[Q^] - F^llQ^]) - ^miQ^))' - (G^fiQ-))') + ^+I/[Q+]^+ 

X < g+, o| p |o, g_ > . (40) 

The non-perturbative correlation function in the presence of the background field A'^^{x) 
is given by 

5 5 



^background QCd[^, J+, V-^ l+i V-]\j+=J- 



=v+=v+=v-=v-=o 



Sr]i{x2) Sr]i{xi) 

= < ijl{x2)i^i{Xi) >background QCD = < in\i:l.{x2)^Pi{Xi)\in >background QCD (41) 

where r, s = +, — are the closed-time indices. The non-perturbative correlation function in 
QCD (without the background field) is given by 

=< ^Ux2)^lixi) >QCD=< tn\^Ux2mixi)\tn >^^° . (42) 
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Using eq. (HOj) we find from eqs. ( HT]) and ( H2|) 

< m|^;(a;2)V^^(Xi)|m >background QCD = 

which gives by using eqs. (fTOjl and ( !T2|) 

< in\ijl{x2WAxi)\in >Lkground QCD = f dX h ■ ^^(x^ + /l^A)T^] 

J — oo 

_ _ r — oo 

X [< m|^;(x2)V5Ka;i)|m >qHd] X [P exp[ig dX h ■ AUx"^ + X)T'']] (44) 

where the closed-time path indices r, s are not summed. All the ^^(x) dependences (respon- 
sible for soft and coUinear divergences) have been factored into the path ordered exponentials 
or Wilson lines. This proves factorization of soft and coUinear divergences in non-equilibrium 
QCD. 



V. CONCLUSIONS 



In order to study hadron production from quark-gluon plasma using the non-equilibrium 
fragmentation function one needs to prove factorization of fragmentation function in non- 
equilibrium QCD. Factorization refers to separation of short distance from long distance 
effects in field theory. The factorization of soft and coUinear divergences in non-equilibrium 
QED is proved in 17|. In this paper we have proved factorization of fragmentation function 
in non-equilibrium QCD by using Schwinger-Keldysh closed-time path integral formalism 
and background field method of QCD in a pure gauge. Our proof is valid in any arbitrary 
gauge fixing parameter a. This may be relevant to study hadron production [9] from quark- 



gluon plasma 



20| at RHIC and LHC. 
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